(N 

o 

(N 

O 
CD 

Q 



CO 



Generalized interface models for transport phenomena: unusual scale effects 
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The effective transport properties of heterogeneous nanoscale materials and structures are affected by several 
geometrical and physical factors. Among them the presence of imperfect interfaces plays a central role 
being often at the origin of the scale effects. To describe real contacts between different phases some classical 
schemes have been introduced in literature, namely the low and the high conducting interface models. Here, we 
introduce a generalized formalism, which is able to take into account the properties of both previous schemes 
and, at the same time, it implements more complex behaviors, already observed in recent investigations. We 
apply our models to the calculation of the effective conductivity in a paradigmatic structure composed of a 
dispersion of particles. In particular we describe the conductivity dependence upon the size of the inclusions 
finding an unusual non-monotone scale effect with a pronounced peak at a given particle size. We introduce 
some intrinsic length scales governing the universal scaling laws. 
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I. INTRODUCTION 

One of the central problems in material science is to 
evaluate the effective electric, magnetic, elastic and ther- 
mal properties governing the physical behavior of hetero- 
geneous materials. 1,2 In recent years, with the progressive 
miniaturization of structures and devices, possible size 
effects have attracted an ever increasing interest. One 
crucial property that usually drives the scale effects in 
structured materials is the complexity of interfaces be- 
tween different phases. Typically, in the macroscopic 
modeling, the interfaces are assumed to be perfect. In 
the context of the electrical conduction it means that the 
potential V and the normal component of the current 
density J are continuous across any interface: 5 ' 4 [V] = 
and [J • n] = 0, where the symbol [/J represents the 
jump of the function / across the interface. This ap- 
proximation turns out to be valid in the case of small 
surface/volume ratio. However, in many real cases of 
technological interest, e.g. nanocomposites, it is impor- 
tant to take into consideration the specific properties of 
the contacts among the constituents. To this aim, two 
effective interface models have been so far introduced for 
describing two extreme situations in a zero thickness for- 
mulation. Moreover, other models are based of an ex- 
plicit interphase of finite thickness and, therefore, they 
typically consider a three-phase heterogeneous material 
composed of the inclusions, the interphase medium and 
the matrix. 

The first zero thickness model is called low conducting 
interface and it is based on the Kapitza resistance, intro- 
duced in the context of the thermal conduction. Accord- 
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ing to this approach [ J ■ rtj = 0, while the potential suf- 
fers a jump proportional to the local flux, [V] = —rJ ■ n, 
where r is the Kapitza-like resistance. The second model, 
called high conducting interface, concerns the case of an 
interphase of very high conductivity with vanishing thick- 
ness. In this situation [VJ = 0, while the normal compo- 
nent of the current density is proportional to the surface 
Laplacian of the potential, [J-n] = gV 2 s V , where g repre- 
sents the interphase conductance. Several investigations 
on heterogeneous materials with low 9-19 or high 12 ' 18 ' 20-24 
conducting interfaces can be found in literature. 

In many cases, the behavior of complex interfaces can- 
not be simply described through the low or high con- 
ducting model. In fact, these schemes account for a sin- 
gle interlayer with an extreme (high or low) value of the 
conductivity, while real interface typically exhibit a com- 
plex or multilayered structure. To overcome this diffi- 
culty we introduce a generalized anisotropic interface for- 
malism, which consider both the normal resistance (simi- 
larly to the Kapitza case) and the tangential conductance 
(as in the high conducting interface model). The term 
anisotropic refers to the fact that the normal resistance 
and the tangential conductance are completely indepen- 
dent, describing a different behavior in the two directions. 
As discussed below, in order to integrate both the normal 
and tangential features, two dual schemes are possible, as 
shown in Fig.l. They exploit the classical T and n elec- 
tric lattice structures. By means of this approach we take 
into account all situations comprised between the low and 
high conducting interface models, which can be seen as 
limiting cases of the present theory. In our schemes both 
the potential and the normal component of the current 
density are discontinuous at the interface. The richness of 
the proposed models allows us to effectively describe the 
behavior of real imperfect/multilaycred/structured inter- 
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FIG. 1. (Color online) Schemes of the dual anisotropic imperfect interfaces (T-model and II-model) between two homogeneous 
media with conductivities a\ and 02. 



faces, which can be found in several heterogeneous mate- 
rials of technological interest and, in particular, for those 
displaying a complex nanoscale structure. 

At first, we have applied the generalized interfaces 
to model a single particle embedded in a different ma- 
trix (inhomogeneity) . One of the most important tech- 
nique used to study this system is based on the Es- 
helby formalism. It has been introduced in the context 
of the isotropic elasticity theory, 25 ' 20 ' generalized to the 
anisotropic elasticity"' and applied to the electric, mag- 
netic or thermal case. 28-30 A relevant universal property 
states that the field induced in cylindrical or spherical 
particles with zero-thickness low or high conducting in- 
terfaces is uniform if the externally applied field is so. 
It is true for both isotropic constituents and anisotropic 
ones, as recently proved. 19 ' 24 In the present work, we 
show that the uniformity property of the internal field 
for spheres and cylinders is preserved also for our gener- 
alized interface models. Moreover, to extend the validity 
of the Eshelby approach, we propose a method which 
is able to determine the field within and around a single 
particle even if the externally applied field is not uniform. 

The previous results have been applied to determine 
the overall conductivity (through an effective medium 
theory ' ) of a dispersion of particles with imperfect 
interfaces. We have verified that, in contrast to per- 
fectly bonded inclusions, the effective properties de- 
pend upon the size of the inhomogeneities. Interest- 
ingly enough, while in the case of low or high conduct- 
ing interfaces the scale effects are described by mono- 
tone scaling laws, the present generalized models show 
non-monotone scale effects with a sizeable peak of the 
transport properties. This point can be considered as 
a specific signature of the complex resistive/conductive 
behavior of the interface. A description of this intrigu- 
ing behavior has been made through different intrin- 



sic length scales governing the universal scaling laws. 
Typical material science problems where these inter- 
face models can be profitably applied are the following: 
tailoring of composites with semiconductor whiskers; 5 5 
thermal optimization of metal/dielectric interfaces' and 
change materials through nanoclusters of stable oxides; 5 1 
analysis of thermal and electric conductivity of carbon- 
based nanostructures; ' size effects understanding in 
SiC/epoxy (or similar) nanocomposites. 38 ' 39 

Throughout all the paper we develop the formalism 
with the terminology of the electrical transport, but all 
results can be applied to the analogous situations of ther- 
mal conduction, antiplanc elasticity, magnetic permeabil- 
ity and electric permittivity as well. 

II. THE DUAL INTERFACE MODELS 

To begin, we introduce a simple lattice network tak- 
ing into account the normal resistors R + and R~ and 
the tangential conductance G (see Fig.l, T- model). This 
structure is able to consider both the anisotropy (along 
the normal and tangential directions) and the different 
behavior of the normal conductivity on the two sides of 
the interface. For the moment we consider a curvilin- 
ear interface between two different materials of a planar 
structure (2D geometry). The generalization to the ar- 
bitrary three-dimensional case will be made straightfor- 
wardly. By a direct application of the Kirchhoff circuit 
laws we obtain two equalities describing the voltage jump 
and current jump across the interface (the definition of 
the relevant quantities is shown in Fig.l) 

v;-Vr = -R+I+-R-I i , (1) 

I+-ir=GR+ (It,- 2lf 

+G {V+, 2V+ + V+J . (2) 
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In the limit of a continuous zero-thickness interface we 
easily obtain from Eqs.(l) and (2) the relations for the 
interface in the form 

= (j-n)", (3) 

I-'"- = (j- H ) + +a^y + - (4) 

The parameters r~ , r + and g are the suitably rescaled 
counterparts of R~ , i? + and G (r~ and r + are mea- 
sured in fim 2 and g in O -1 ). In previous expressions, 
the partial derivatives are performed with respect to the 
variable s, which represents the curvilinear abscissa along 
the arbitrarily curved interface on the plane. As usual, in 
the three-dimensional case the operator d 2 /ds 2 must be 
substituted with the surface Laplacian V|, which is in- 
troduced and discussed in Appendix A. We can observe 
that the present approach reproduces the low conduct- 
ing interface model if g = (with a Kapitza resistance 
r = r~ + r + ) and the high conducting interface model if 
r~ = r + = 0. 

The low conducting model is characterized by a se- 
quence of normal resistances R = R + + R" (T-modcl 
with G = 0). If we consider As as the step along the 
curvilinear abscissa s and Az as the step along the direc- 
tion perpendicular to the plane represented in Fig.l we 
have R = ^stsi where 

a± is the normal conductivity 
of the interphase of thickness Ah. The Kapitza resis- 
tance is therefore given by r = RAS where AS" = AsAz 
is the area element associated to J ■ n; we finally obtain 
r = lim / \f l ~>o,cr ± ~¥0 ^- 12 Similarly, the high conducting 
model is characterized by a series of tangential conduc- 
tances G (T-modcl with R~ = R + = 0). It is simple 
to observe that G = a\\ A, ^ z where o\\ is the tangential 
conductivity of the interphase. The specific conductivity 
is therefore given by g — G^r (where AS" = AsAz) and 
we obtain the result g = liniAh^o,<T||^oo cr\\Ah. 12 So, we 
have a direct link between the interphase properties (Ah, 
ax, crii ) and the models parameters (r, g) for the high and 
low conducting interfaces. Interestingly enough, we ob- 
serve that when we consider an anisotropic single layer in- 
terphase (which is uniaxial or transversely isotropic with 
normal conductivity a± and tangential conductivity cry), 
the only component aj_ is relevant for the low conducting 
model and the only component o~\\ is relevant for the high 
conducting interface. 

Of course, if both relations g = and r~ = r + = 
are satisfied in the T-model, then the ideal inter- 
face is simply obtained. It is not difficult to prove 
that this model is completely equivalent to a series 
of three different ideal sheets (multi-layered interface) 
A, B and C: an external low conducting phase with 
Kapitza resistance ? ,+ = lim Afe(A:) _^ A ^ , a 

halfway high conducting phase with specific conduc- 
tance g = lim^^( B )_^g o-tB) -too A/i^" 8 ) and, finally, an 
internal low conducting phase with Kapitza resistance 



r = lim A/.(O)->0,<rf ^0 l^T - Tlie three la y erS are 

characterized by thickness AhS A \ Ah^ B \ Ah^ (with 
Ah = AhW + Ah^ + Ah^) and conductivities a {A \ 
o~^ B \ 0j_ . So, we have built an example of interpreta- 
tion of the model parameters with a concrete physical 
multilayered structure. 

A dual model can be introduced by considering the 
second structure depicted in Fig.l (II-modcl). A proce- 
dure similar to the previous one leads to the following 
interface equations 

M = ~r (/• n) + + rg + -^V+, (5) 

where the parameters r, g + and g~ are the suit- 
ably rescaled counterparts of R, G + and G~ , appear- 
ing in Fig.l, right. As before, the operator d 2 /ds 2 
must be substituted with the surface Laplacian V| 
for the 3D case. We can prove that also the II- 
modcl is exactly equivalent to a series of three differ- 
ent ideal sheets: an external high conducting phase 
with conductance q + = lim.,,,, _ <a> a 

^ Aru- 4 '— >0,(7|| — >oo || ' 

a halfway low conducting phase with Kapitza resis- 
tance r = lim A , (B j „ a <.B) ^tbt~ an d, finally, an in- 

' _L 0j_ 

tcrnal high conducting phase with conductivity g~ = 
lim.,™ _ (o ai ^ Ah^ c \ As before, the layers are 

Ah^J— t0,a:, '-too \\ ' J 

characterized by thickness Ah^- A \ Ah^ B \ Ah^ and con- 
ductivities o~^ A \ o~^\ o"jj • It is important to remark 
that the interpretation of the model through three ad- 
jacent layers (for both the T and II structures) it is not 
restrictive; in fact, the proposed schemes can be also used 
to effectively represent different imperfect interfaces with 
all parameters fitted in order to mimic their correct be- 
havior. We also underline that some more complete mod- 
els have been proposed in literature (see, e.g., the recent 
Gu and He interface,' which also degenerates to the high 
or low conducting models and it is able to take into ac- 
count all the coupling among the electric, magnetic and 
clastic fields); here, we have proposed our schemes with 
the idea to find a compromise between the complexity 
and the possibility to analytically solve the problem for 
paradigmatic composite structures. 

The proposed models are dual from both the geomet- 
rical point of view, as shown by the T and II lattice 
structures, and the physical point of view, as discussed 
below through the results for the composite materials. 



III. SINGLE PARTICLE BEHAVIOR 

We consider now a single circular (in 2D) or spheri- 
cal (in 3D) particle with conductivity embedded into 
a matrix with conductivity o\ (see Fig. 2): we suppose 
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FIG. 2. (Color online) Scheme of a single circular (a) or spher- 
ical (b) particle with conductivity 02 embedded into a matrix 
with conductivity 01. The interface between the two phases 
is described by either Eqs.(3) and (4) or Eqs.(5) and (6). 



that the interface between the constituents is described 
by Eqs.(3) and (4) (T- model) and we determine the effect 
of an arbitrary externally applied field. Since we are deal- 
ing with two isotropic phases, in order to solve the prob- 
lem we can directly apply the original idea of Maxwell, 111 
which is based on the following steps. Firstly, we ob- 
serve that the electrical potential must be an harmonic 
function both inside and outside the particle: therefore, 
it can be straightforwardly expanded in trigonometric 
series (in 2D) and in series of spherical harmonics (in 
3D). 41 Secondly, we can substitute such expansions in 
the interface conditions, by obtaining a set of equations 
for the unknown coefficients, completely describing the 
potential both inside and outside the inhomogeneity. To 
accomplish this last step we must determine the surface 
Laplacian of the series expansions: to do this we remem- 
ber that the trigonometric functions and the spherical 
harmonics are eigenfunctions of the V| operator with 
certain eigenvalues described in Appendix A. The com- 
plete procedure, which is valid for any externally applied 
field, is described in Appendix B for the 2D case and in 
Appendix C for the 3D case. Here, we are interested in 
the particular case with an uniform applied electric field 
Eq, corresponding to a potential Vq = —pcostfEn (see 
Fig. 2). The perturbation induced by the inhomogeneity 
with imperfect contact has been eventually found as 



for p < R V = -pcostfEo 



da 1 



(7) 



for p > i?=> V = -pcostiEn [ 1 + ^j- ] , (8) 



R*_B 
p d C 



where d = 2 for the circle, d = 3 for the sphere and the 
parameters B and C are defined as follows 



B = o\ — 02 + - 

-(d-l)i- 
V ' R 



R 



-01CT2 



1 — r T — 
R 



01 

1 + r ~R. 



C = (d- l)o-i + <J2 + (d-l) 
9 



R 



-<J\<J2 



+(d-l)^[l + (d-l)r 



11 
R 



(9) 



.(10) 



We can observe that the electric quantities both inside 
and outside the particle, in contrast to the case with per- 
fect interfaces, depend on R. So, the previous result can 
be used for analysing the scale effects induced by the 
imperfect contact. From Eq.(7) it is easy to identify the 
induced internal field as Ei nt / Eq = d<j\jC. A first scaling 
law for R — > 00 can be obtained by introducing the clas- 
sical Lorentz field for a particle with a perfect interface 



E u 



Ej, 



r + =r ,- = o,o=o; we can easily prove that 



Ej. 



(d-l)cri £-+£+ + £ 



Eio 



(d — l)a\ + (72 



R 







1 

R? 



(11) 



where we have introduced the following intrinsic length 
scales 

r =a 2 r~, £+ =a 2 r + , C = (12) 

0"! 

which automatically emerge from the analysis and com- 
pletely control all the scaling laws. Eq.(ll) means that 
the internal field approaches the Lorentz field for large 
radius of the particle (R » l~ + £ + + £), i.e. the ef- 
fects of the contact imperfection are vanishingly small 
for R — > 00. 

We discuss now the scaling laws obtained for R —> 0. 
A long but straightforward analysis leads to 



E,, 



doo 



En 



Eint 



En 



3=0 



Ei 



En 



(d- 1) 2 cti £-£+C 
d R 

(d -!)£-+ £+ + 

d R 



r+ —r~ — 



{d-l)C 



OR' 



+ 0(R 4 ), (13) 
O (i? 2 ) , (14) 
(15) 



In any case the internal field converges to zero for very 
small particles. It is interesting to observe that the in- 
ternal field for the T-model follows a scaling law with a 
power of three, while the low and high conductivity mod- 
els follow a law with a scaling exponent equal to one. It 
can be seen in Fig. 3 where log\n{Ei n t / Eq) is represented 
versus log\nR. The blue curves (with squares) and the 
black ones (without symbols) describe the generic inter- 
face (g 0, r + 7^ 0, r~ 7^ 0) and show a slope +3 for 
small R, which is in agreement with Eq.(13). On the 
other hand, green curves (with triangles) and red ones 
(with circles) correspond to the high and the low con- 
ductivity interface, respectively: they all exhibit a slope 
+1 for small R as predicted by Eqs.(14) and (15). We 
also note that all curves in Fig. 3 converge to the Lorentz 
field for R — > 00, as described by Eq.(ll). 

Now, we take into consideration the LT-model described 
by Eqs.(5) and (6). The perturbation to the electric po- 
tential generated by the inhomogeneity is described again 
by Eqs.(7) and (8) but with new coefficients B and C 
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(15) become as follows 

Eint da 2 



Eq 



En 



Eint 



(d - 1)2(71 ec+c- 



O (i? 4 ) , (20) 



R 



(d- 1) t 
d R 



r=0 



(d -!)£+ + £ 



0{R 2 ), (21) 
O (R 2 ) . (22) 



As before, the internal field converges to zero for very 
small particles (with different scaling exponents, as above 
described). 



IV. EFFECTIVE CONDUCTIVITY OF DISPERSIONS 



FIG. 3. (Color online) Plot of log w {E int / E ) ver- 
sus logioR for the T-model. Green curves with trian- 
gles: high conductivity model with a varying g in SI = 
{0.001,0.016,0.25,4,64,1000}. Red curves with circles: low 
conductivity model with a varying r + — r~ in fi. Blue curves 
with squares: general model with r + — r~ = 1 and g vary- 
ing in fi. Black curves without symbols: general model with 
g = 1 and r + = r~ varying in tt. Everywhere, the dashed 
lines correspond to the values < 1 of the varying quantity. 
Parameters adopted in a.u.: a\ = 1, <T2 = 5, d = 3, c = 0.3. 



given below 

B = o\ — (72 4 

(d-1) 
R 



r 2 9 9 +r 



R 



,+ 



1 + r 



11 
R 



}, (16) 



C = (d - 1)^7! + O-a + (d - 1)^(7^2 + (d - l) 2 ^^ 

R R 6 



R 



l + (d-l)r^' 



1 , CT2 
1 + r i? 



(17) 



With regards to the scaling law for R — > oo, it is possible 
to prove that Eq.(ll) must be substituted with 



Eint 
Elor 



- 1 = - 



(d - l)(7l + (7 2 



R 



(18) 



To analyse the effects of imperfect interfaces on a com- 
posite material we consider a dispersion of cylindrical or 
spherical particles of conductivity 02 in a matrix with 
conductivity <7i. When the interfaces are described by 
Eqs.(3) and (4) (T-model) we can generalize the Maxwell 
approach 1 " or, equivalently, the Mori-Tanaka scheme 42 
by obtaining the following effective conductivity for a 
composite with a volume fraction c of the dispersed par- 
ticles 



J eff 



f 



cdB 



(23) 



{l-c)C + c[C- (d-l)B] 



where B and C are given in Eqs.(9) and (10). Detailed de- 
scriptions of the homogenization procedures can be found 
elsewhere. 19 ' 24,31 ' 32 For interfaces described by the low 
conductivity model we obtain ai ow = cr e // |g=o, which 
is in perfect agreement with recent investigations; 18 ' 19 
on the other hand, when the high conductivity model is 
accounted for we have Ohigh — v e ff |r+=r-=o> wn ich cor- 
responds to some known results. 18 ' 24 Moreover, when we 
consider a perfect contact between the constituents we 
obtain the celebrated Maxwell formula 40 



1 



1 



dc((7l - (7 2 ) 



(24) 



(1 



[(d - l)(7l + (7 2 ] + cdcr 2 



The first important scaling law concerns the situation 
with a large radius of the particles: in this above 
said, the size effects disappear and the effective conduc- 
tivity converges to the Maxwell one as follows 



where we have introduced the dual intrinsic length scales 



- 1 = 



cd 2 aj n 
Q R 



O 



1 



(25) 



e = (7 2 r, C+ = y —, C- = y —. (19) 

(7l (7l 

As expected, also in this case the internal field approaches 
the Lorentz field for large radius of the particle (R 3> 
I + £+ + C~). On the other hand, for R -> 0, Eqs.(13)- 



where we have defined 

n = (d-l)£-^(£+ +r), (26) 
0-1 

= [(d + c-l)(7i + (l- C )c7 2 ] 

x [(d-l)(l-c)(7i + (cd-c+l)<7 2 ]. (27) 
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The parameter % represents the overall length scale of 
this process and it is a linear combination of the terms 
defined in Eq.(12). This result can be simply compared 
with recent achievements 1 ^" concerning the cases with low 
and high conductivity interfaces. Indeed, we find a per- 
fect agreement if l + = l~ = or C = 0. 

Other interesting scaling laws can be found for R —> 0. 
To analyse this case we define the conductivity 00, which 
represents a Maxwell dispersion with a 2 (dispersion 
of voids), and the conductivity CToo, which characterizes 
a Maxwell dispersion with 02 — > 00 (dispersion of super- 
conducting particles): 



01 



(1-cXd-l) , 
d+c-1 ' 



"00 

01 



1 



cd 



1 



(28) 



First of all we observe that if r + =/= we have cr e ff —t&a 
with the scaling law 



a eff 

00 



1 



cd 2 a 2 



R 



(d-l)(l-c)(d-l + c)o-i £+ 



Similarly, if r ^ with r + = 
u e ff — >0o with the scaling law 



and g 



+ O (R 2 ) . 

(29) 

we obtain 



a eff 

00 



1 



Cd 2 (72 



R 



(d-l)(l-c)(d-l + c)o-i t 



O (R 



(30) 



So, for the general T-modcl and for the low conducting 
interface we have cr e ft — ><tq (when R 0) with a scaling 
exponent equals to one. On the other hand, for g 



and r + = we prove the convergence <7 e // 
scaling law 



with a 



T eff 



- 1 = - 



cd 2 



R 



(d-l)(l-c)(cd-c+l) C 



O (R 



(31) 



It means that the high conductivity model leads to 
(Jeff -^-cfoo for R ->• 0. 

This complex scenario is summarized in Fig. 4 where 
e // is shown versus logioR. Even at constant volume 
fraction c, significant size effects on the effective conduc- 
tivity are evident for a variable radius R. In Fig. 4 (top) 
we have reported the results for 02/ 1 <J\ = 2 and in Fig. 4 
(bottom) for (72/01 = 0.5. In both cases we have shown 
the neutrality axis at which the effective conductivity 
cr e ff equals the matrix conductivity 01, making the in- 
clusions effectively hidden. 12 ' 4 '' We can observe that a e ff 
is a monotonically decreasing function of R (from Uoo to 
&max) for the high conductivity model (green curves with 
triangles), while it is a monotonically increasing func- 
tion of R (from oq to a max ) for the low conductivity 
model (red lines with circles). So, the neutrality con- 
dition (B = 0) can be satisfied by the low conductiv- 
ity model for 02 > 01 (02 — 01 = Ttj-iaifR, sec Fig. 4, 
top) and by the high conductivity model for 02 < 01 
(ax — 02 = g(d— 1)/-R, see Fig. 4, bottom). On the other 
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FIG. 4. (Color online) Plot of er e // versus logioR for the 
T- model. We adopted the following parameters (in a.u.): 
01 = 1, <J% = 2 (top) and 01 = 1, 02 = 1/2 (bottom). Every- 
where we used d — 2, c = 0.3. Green curves with triangles: 
high conductivity model (r + = r~ = 0) with a varying g 

in n = |l0" 3+2(j " 1)/3 ,j = l...lo}. Red curves with circles: 

low conductivity model (g = 0) with a varying r + = r~ in SI. 
Blue curves with squares: T-model with r + = r~ — 1 and p 
varying in SI. Black curves without symbols: T-model with 
g = 1 and r + = r~ varying in SI. Everywhere, the dashed 
lines correspond to values < 1 of the varying quantity. 



hand, the blue and black lines concern the case of the gen- 
eral T-modcl and they exhibit a non monotone behavior 
starting from 00 and arriving at a max . It is interesting to 
note that, with the general T-model, it is possible to sat- 
isfy the neutrality condition for both the cases 02 > 01 
and 02 < 01 • The condition leading to neutrality in this 
case (T-modcl) is 



01 - 02 



0102 



(32) 



(d-l)i[l-r+t] + 
which is represented in Fig. 4 by the intersections of blue 



7 




1.5 



a eff 


log w R 1 


Coo 


\\\Yi\\V 




\\\ \\ \\ \ 
■ \ ' + \ • \ 

\ V- V- '>'', 
\ ' ; \^ V> 

neutrality, ^ \ 


\ \\\\ \ \ \ \ 

A \ ft. \ \ \ \ 
^YxVb \ \ \ \ 


_i\ V ^y^^'yy' 17 max 
,-):,i:Mi-f,-t'd4^Vf>^ 



'-5 



log 1Q R 



FIG. 5. (Color online) Plot of (Jeff versus logioR for the 
Il-model. We adopted the following parameters (in a.u.): 
<7i = 1, (J2 = 2 (top) and oi = 1, 02 = 1/2 (bottom). Every- 
where we used d = 2, c = 0.3. Green curves with triangles: 
high conductivity model (r = 0) with a varying g~ — g + 

in Q. = |l0 _3+2(:,_1)/3 ,j = 1...10|. Red curves with circles: 

low conductivity model (g + — g~ — 0) with a varying r in fi. 
Blue curves with squares: Il-model with r = 1 and g + = g~ 
varying in Q. Black curves without symbols: Il-model with 
g + = g~ — 1 and r varying in 17. Everywhere, the dashed 
lines correspond to values < 1 of the varying quantity. 



and black curves with the neutrality axis. 

As for the dual Il-model, we can affirm that the gener- 
alized Maxwell theory given in Eq.(23) is still valid but 
the coefficients B and C must be taken from Eqs.(16) and 
(17), respectively. For a large radius of the particle we 
have the scaling law identical to Eq.(25) where Q is given 
by Eq.(27) and H by the following expression 



n = (d -!)(£+ + C~) 



01 



(33) 



It represents the overall length scale of the Il-model and 
it is indeed a linear combination of the terms defined in 



Eq.(19). We also report the scaling laws for R 0. If 
g + 7^ we have that e // — >(Too with the scaling law 



a eff 



- 1 = - 



cd 2 



R 



(d- 1)(1- c){cd- c+1) £+ 



OR 2 



(34) 



Similarly, if ^ with r = and g + = we obtain 
<7 e / / —>cr 00 with the scaling law 



J eff 



- 1 = - 



cd 2 



i? 



(d-l)(l-c)(cd-c+l) £- 



O (R 2 ) . 
(35) 



So, for the general il-model and for the high conducting 
interface we have tr e / / — >croo (when R — > 0) with a scaling 
exponent equals to one. On the other hand, for r =/= 
and g + = we prove the convergence o^// — >cro with a 
scaling law 



0Q 



1 



cd 2 cr 2 



1? 



(d-l)(l-c)(d-l + c)cri I 



0(R 2 



(36) 



It means that, as expected, the low conductivity model 
leads to cr e /f — >-0o f° r -R — > 0. 

In Fig. 5 the results for the Il-model are shown: the 
effective conductivity is represented versus the radius R 
of the particles. In Fig. 5 (top) we have the case with 
02/01 = 2 and in Fig. 5 (bottom) we show the results 
for (72/01 = 0.5. All previous scaling laws are confirmed 
and clearly indicated. By drawing a comparison between 
Fig. 4 and Fig. 5 we can point out the dual character of 
the proposed models: the T-model behaves similarly to 
the low conducting interface with regards to the limiting 
cases R — > and R — > 00 but it shows a specific ad- 
ditional upwards peak describing the competition of the 
scale effects with the presence of the tangential conduc- 
tances g. Conversely, the Il-model behaves similarly to 
the high conducting interface with regards to the limiting 
cases R — > and R — > 00 but it shows a specific addi- 
tional downwards peak describing the competition of the 
scale effects with the presence of the normal resistances 
r. As before, also in the case of the Il-model, we can 
satisfy the neutrality condition for both the contrast sit- 
uations 02/01 > 1 and (72 /oi < 1. The condition leading 
to neutrality in this case (Il-model) is 



(7i - cr 2 



g + +g~ 

R 



(d-1) 



(d-1) 



2 g + +g 
R 3 



d Kr(9-<ri - .9+0-2) 



R 



,(37) 



and it is satisfied in Fig. 5 at the intersection points be- 
tween the blue or blacks curves and the neutrality axis. 

By means of this analysis we can assert that the T 
and n models exhibit an interesting complex behavior 
which is able to reproduce many properties of real inter- 
faces appearing in different nano-systems. As an example 
we can compare our results with those recently obtained 
for a dispersion of SiC particles (with radius between 5 
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and 15A) in a polymeric (epoxy) matrix. 39 By means of 
a multiscalc combination of the non-equilibrium molec- 
ular dynamics and a micromcchanics bridging model, 
the thermal conductivity has been studied in terms of 
the particles radius. The result is in perfect qualitative 
agreement with our T-model and a maximum value of 
the conductivity was obtained for a given radius. To 
obtain such a result the Kapitza resistance and a spe- 
cific interphase describing the bonding of the polymers to 
the monocrystallinc SiC particles have been considered.' 9 
Our T-model is able to describe the overall response of 
the structured/multilayercd interface through the simple 
conditions given in Eqs.(3) and (4) imposing the jumps 
of the physical fields over the zero-thickness interface. 
Therefore, the proposed models perfectly implement the 
multiscalc paradigm by introducing the effective proper- 
ties of a given interface behavior. 

We remark that in this Section we have used the gener- 
alization of the Maxwell approach 40 or the Mori-Tanaka 
scheme 42 in order to obtain simple results and to directly 
analyse the scale effects induced by the imperfect inter- 
faces. Nevertheless, the closed form results discussed in 
Section III for the single particle response can be easily 
exploited to implement other homogenization techniques 
such as the differential method 31 ' 44 ' 45 , the self consistent 
scheme 46-48 , the generalized-self-consistent model 49 and 
the strong-property- fluctuation theory. We also remark 
that the analysis of the imperfect interfaces is an impor- 
tant topic also in the field of micromechanics (elastic- 
ity of composites) where several theoretical models have 
been proposed 51-5 ! and intriguing scale effects have been 
observed. 53-55 



V. SUMMARY AND CONCLUSIONS 

In this paper we have taken into consideration the pos- 
sible scale effects induced by imperfect interfaces between 
the constituents of an heterogeneous system. To this aim 
we introduced two generalised schemes, namely the T and 
II structures, which can be seen as natural combinations 
of the so-called low and high conducting interface mod- 
els. One important property discussed concerns the uni- 
formity of the physical fields in circular or spherical parti- 
cles with T or II imperfect interfaces. This point extends 
well known theorems proving the uniformity in different 
conditions and opens the possibility to study the behav- 
ior of new interfaces in anisotropic, elliptic and ellipsoidal 
particles, which are standard problems in the theory of 
inhomogeneities. The results for a single inclusion were 
applied to the analysis of the effective properties of dis- 
persions. In particular we studied the scale effects and 
we found interesting behaviors, which generalize those 
observed with low and high conductivity interfaces. We 
indeed observed a specific peak of the effective conductiv- 
ity in correspondence to a critical radius of the dispersed 
particles: it corresponds to the competition between the 
tendency to attain the Maxwell conductivity limit for a 



large radius and the conduction properties of the inter- 
face, which tend to increase or decrease the overall con- 
ductivity, depending on the specific parameters. This is 
exactly the trend observed in recent analysis of imperfect 
interfaces in nanocomposites (hard particles in polymeric 
matrix or similar mixtures). To conclude, we have anal- 
ysed the neutrality properties of the T and II models: 
contrarily to the low and high conducting interface, we 
have proved that it is possible to satisfy the neutrality 
condition for any contrast <Jij o\ between the conductiv- 
ities of the involved phases. So, Eqs.(32) and (37) are the 
updated versions of the neutrality criteria, representing 
the generalizations of some findings, published in recent 
literature. We remark that all the achievements of the 
present paper can be also used in dynamic regime if we 
consider a wavelength A of the propagating waves that is 
much larger than the radius R of the particles. In this 
case we are working in the so-called quasi-static regime 
and any inhomogeneity feels a nearly static applied field. 



Appendix A: The surface Laplacian 



The surface Laplacian operator is defined as 



1 d 

y/gdati 



2L 

dcti 



(Al) 



where gij are the components of the metric tensor (the 
first fundamental form) of the Ricmannian manifold (the 
surface) r = f(ai, aa)- 56 It means that g^ = • -J^- 

and the dual components g lJ are obtained by inverting 
the matrix g^ . The quantity g is the determinant of g^ . 
Typically, in differential geometry of two-dimensional 
surfaces we adopt the symbols gn = E, g±2 = gi\ = F 
and 1722 = G; so, for an orthogonal system of coordinate 
lines F = and Eq.(Al) reduces to 




Gdf_ 
E dai 



d 



Edf_ 

G dao 



(A2) 



For a planar circle r = (i?cos?9, i?sin$) we simply 
have 



V 2 f = ^l = J_^l 

sJ ds 2 R 2 dd 2 ' 

and the following property is evident 

1 



R 2 



n 2 e in ». 



(A3) 



(A4) 



It means that the trigonometric functions cos 71$ and 
sin rvd are eigenfunctions of the Laplacian operator with 
eigenvalues — jpn 2 - 

For a spherical surface r = (iJcos^sini?, i?sin(/?shn9, 
Rcos-d) it is possible to obtain 



1 



1 d 



R 2 smddd 



1 d 2 f~) 



sin 1? dip 2 
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and we can prove that 

V 2 s Y nm (ti,ip) = -ifi(n + l)F„ m (i), V ). (A6) 

It means that the spherical harmonics Ynmid, if) are 
cigenfunctions of the surface Laplacian operator with 
eigenvalues — jpn(n + l). 41 They are defined (for n > 0, 
— n < m < n) as ' 



where P n n (£.) are the associated Legcndrc 
polynomials 57 ' 58 

m 1 jn+m 

p: i (o = (i - ey ^_de^ ( e - 1)" .(as) 



Appendix B: Two-dimensional geometry: the circle 

We suppose to consider a circular inhomogeneity of 
radius R (conductivity a 2) in the plane (x,y) with con- 
ductivity <j\. We consider an arbitrary applied (or pre- 
existing) potential Vq(x, y) and we search for the pertur- 
bation induced by the inhomogeneity. Since the electric 
potential must be harmonic both inside and outside the 
interface, we have 

+00 

V = V + P n ( A n cosnti + B n sinn??) ,p < R, (Bl) 

n=0 

+ 00 

V = Vb + ^ P~ n (An cos nd + B n sin nv) ,p> R, 

where (p, •&) are the standard polar coordinates. The 
potential Vq and its derivatives 4^ can be expanded in 
Fourier series for p = R 

+00 

V (R,-d) = J2(C n cos nd + D n sin 711?) , (B2) 

n=0 

dV + °° 

(R, ■&) = (F n cos nQ + G n sin nff) . 

P n=0 

By substituting Eq.(Bl) in the anisotropic interface 
model (Eqs.(3) and (4)) and by using Eqs.(A4) and (B2), 
we obtain a set of equations for A n and A n 



R- n A n - R n A n = r+<n [F n - nR~ n - l A r 
+r-a 2 (F n +nR n - 1 A n ) , (B3) 
a 2 (F n + nR n - x A n ) - a x (f„ - nR-^A.^ 

= gr + a x R- 2 n 2 (f„ - nR-^A^j 

-gR- 2 n 2 (C n + R- n A^j , (B4) 



and a similar one for the unknowns B n and B ni not 
reported here for brevity. These systems can be eas- 
ily solved obtaining the electrical potential in the whole 
plane. In the particular case of a uniform applied field 
Vq = —xEq = —pcos{)Eo only the coefficients A\ and 
A\ are different from zero and we obtain Eqs.(7) and (8) 
for d = 2. A similar procedure (not reported here for 
brevity) can be followed to analyse the properties of the 
Il-model described by Eqs.(5) and (6). 



Appendix C: Three-dimensional geometry: the sphere 

We consider now a spherical inhomogeneity of radius 
R (conductivity er 2 ) m a matrix with conductivity a±. 
As before, we assume an arbitrary applied potential 
Vo(x, y, z) and we study the effects of the embedded par- 
ticle. The final electric potential can be expanded as 
follows 

+ 00 +n 

V = V + Y J B nmP n Y nm {d^),p<R, (CI) 

n— m——n 
+ 00 +n 

V = V + J2 C nm p- n - 1 Y nm (-3,ip),p>R, 

n—0 m=—n 

where we have introduced the spherical coordinates 
(p,d,ip). The potential Vb and its derivatives can 
be expanded in a series of spherical harmonics for p = R 

+ OC +71 

V (R,&,<p) = J2 £ PnmYnm($,V), (C2) 



n—0 m— — n 

+ 00 +n 



dp 



Y nm ($,<p). 



n—0 m— — n 



By substituting Eq.(Cl) in the anisotropic interface 
model (Eqs.(3) and (4)) and by using Eqs.(A6) and (C2), 
we obtain a set of equations for B nm and C nm 



R~ 



= r+a l [a nm - (n + l)Rr n ~ 2 C nm ] 

+r~a 2 [anm +nR n ~ 1 B nm ] , 

0"2 [ a nm + nR n ~ 1 B nm ] 

-0-1 [a nm -(n + l)R- n - 2 C nm ] 
= gr + a x R- 2 n{n + 1) [a nm - (n + l)R,^ 2 C nm ] 



(C3) 



-gR- 2 n{n + 1) [p nm + iT^C,, 



(C4) 



It is now possible to find B nm and C nm obtaining the 
electrical potential in the whole space. In the particular 
case of a uniform applied field Vq = — zEq = —pcosfiEo 
only the coefficients B\q and C10 are different from zero 
and we obtain Eqs.(7) and (8) for d = 3. We remark that 
a similar procedure can be followed for studying the dual 
interface described by Eqs.(5) and (6). 
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